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ABSTRACT 



A one dimensional motion of the Bethe-Johnson gas is studied in a context of 
Landau's hydrodynamical model of a nucleus-nucleon collision. The expressions for 
the entropy change, representing a generalization of the previously known results, 
are found. It is shown that these expressions strongly depend on an equation of 
state for the baryonic matter. 

1 INTRODUCTION 

The application of special relativistic gasdynamics to some problems in par- 
ticle physics goes back to the pioneering work by L. Landau P who obtained 
the asymptotic solution describing a one- dimensional expansion of a finite 
slab of matter modelled by the ultra-relativistic gas with the adiabatic index 
7 = 4/3. Later LKhalatnikov 2] obtained the exact solution of the same 
problem. G.Milekhin [3j extended this solution to an ultra-relativistic gas 
with an arbitrary adiabatic index. S.Belenkiy and G.Milekhin |lj applied a 
model of the ultra-relativistic shock wave to a problem of an entropy increase 
in a high energy nucleus-nucleon collision. A significant input into the prob- 
lems associated with a hydrodynamical model of nucleus-nucleon collisions 
was made by G.Baym [5^ who solved the Landau problem [Ij for the cylin- 
drical geometry. Nonetheless G.Baym and collaborators indicated that there 
is a need to use a more realistic equation of state than the one with F = 4/3. 
Since there is still some uncertainty about an exact form of the equation of 
state for the baryonic matter at very high densities, one can consider a wide 
range of possible candidates for such an equation. 

One of these equations was proposed by H.Bethe and M.Johnson [HI (which 
in what follows we call BJ equation), and later it was used to illustrate an 
idea of a possible baryon-quark transition [7j. This equation was also used 
for a speculation on a potential influence of such a transition on a collec- 
tive flow in a high energy ion collisions [HI. Therefore, it seems reasonable 
that a study of one-dimensional motions of Bethe-Johnson gas( including 
the original Landau problem) might provide a new insight into high energy 
phenomena within a model based on relativistic gasdynamics. 
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2 SHOCK ADIABATICS OF BJ GAS 



We begin by considering a relativistic shock. In what follows we adopt the 
units with the speed of light taken to be c = 1. The energy, momentum, 
and particle conservation laws across the shock (in the shock frame)yield the 
following shock adiabatic jH] 

f^l-f^2 + liPl-P2)iVl + V2) = (1) 

where indices 1 and 2 denote the areas in front and the back of the shock 
respectively, /i = W = w/n, w = e + p is the enthalpy per unit of 

volume per unit of baryon rest mass, n is the particle density, p is the pressure 
per unit of baryon rest mass, e is the energy density per unit of baryon rest 
mass, 

V = W/n = wjn^ (2) 

is the generalized "specific volume" ^U]. In these variables the relativistic 
shock adiabatic looks exactly as the Hugoniot adiabatic of a newtonian 
gas. 



2.1 SOME RELATIONS FOR A RELATIVISTIC SHOCK 

We provide a series of relations which follow from the conservation laws across 
the shock. From the continuity of particle flux density 

and momentum flux density 

(-)lf/l +Pl = (-)2f/2'+P2 

n n 



(where we use the definition W = w/n) we obtain 

P2 -Pi 



J' = (3) 



indicating that either 



«)P2 > Pl, Vi > V2 

b)p2 < Pl, Vi<V2 (4) 
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in exact analogy to the newtonian case. It has turned out that only case a) 
can actually occur. 



This is easily shown by considering a relativistic weak shock wave, in which 
the discontinuity in every quantity is small as compared to the quantity. Let 
us expand both sides of the shock adiabatic (P) in powers of small differences 
of entropy(per unit of volume) As = S2 — Si and pressure Ap = p2 — Vi- 
Following the analogous procedure in newtonian gasdynamics, the expansion 
in Ap is to be carried as far as the third order. Thus we obtain: 

/il - /i2 = - Wl = + (^)iA p + 

where according to the thermodynamic identities, 

,dW, _ , , , , , ,dW, 

(--—)„ = T [the absolute temperature), (— — = V 
OS op 

Therefore (0) yields: 

dV 

/xi - /i2 = 2W^T^/\s + 2W^V^/\p + [V^ + ),](Ap)2 + 

Opi 

dV d'^V 

The generalized "specific volume" V = WV needs to be expanded only 
in terms of Ap, since the terms with As contribute to the fourth order of 
smallness in Ap. Thus we get after some algebra: 

V,-V, = W2V2-W,Vi = [V,' + W,&s]{Ap)\Ap) + l[3V&s]Ap' (7) 

Opi 2 Opi 

where once again we use the thermodynamic identity V = {dW/dp)s. Sub- 
stituting this expansion together with (jH)) in ([T]). we obtain 
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which once gain represents the exact analogy of the respective relation in 
newtonian gasdynamics (e.g., Ref. P). In all cases that have been studied 
the second derivative 



> 



dpi 

Therefore by the law of increase of entropy S2 > Si and according to ^ 
P2 > Pi, which implies (EqE)) Vi > V2. 

It was also shown [H] that a weak relativistic shock waves propagates with a 
speed of sound = {dp/de)s- In fact, for small p2 — pi and V2 — Vi relation 
(0) can be written, in the first approximation, as 

2 r9p. 



J 



In the same approximation 



[/2 = t/2 = u^ = jV^ = -V\^)s (9) 

dV 

On the other hand, 

dp_ ^ dp 1 
^dV^' ^de^'[{dV/de)s 
where by using the identity 

we find 

dV _ _ 1 dn _ _V _ _ 1 
oe oe w nw 

and 

/OT^/ON fdvoV"^, j,2t9(e + p), ,dV. 
{dV/de)s = = V^l ]s + 2Vw{—)s = 

As a result Eq.® yields 



(10) 
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Figure 1: Pressure behind the shock p2 as a function of V2/V1 for 1/ = 1.3 and a 
cold gas {pi = 0) in front of the shock 



2.2 SHOCK WAVE IN BJ GAS 

We consider a shock wave propagating in a specific gas with the following 
equation of state (Bethe- Johnson gas) 

p= (zy- l)(e-n) (11) 

where u is an arbitrary constant playing the role of the adiabatic index. If 
the gas in front of the shock is cold, that is pi = , equation (fTT|l yields for 
that region: 

ei = ni (12) 

which implies /xi = 1 and hence Vi = Vi, that is the conventional specific 
volume. 

By substituting Eqs. (fTT|) . (fT^ in the shock adiabatic (P) we obtain the 
following equation with respect to p2 

P2[P2V,{^^ - ir - ^^)(^ + 1)^ + 1] = (13) 

It has two solutions , one corresponding to a trivial case of a continuous flow 
with P2 = Pi = ( weak shock in this case), and the second one representing 
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the strong shock wave: 



\V2 - Vi{v - i)Y 

A graph of p2^i vs. V2/V1 for a specific value of the "adiabatic index" z/ = 1.3 
is shown in Fig.l. 

From equation ()14j) follows that it has two branches ( as is clearly seen 
in the Fig.l), one for V2/V1 < u — 1 and another one for V2/V1 > u — 1. 
We cannot determine the physical branch on the basis of the adiabatic only 
, since for both branches the conditions p2 > Pi,V2 < V\ [which follow 
from the law of entropy increase (jH)) and Eq. ljllll ] are satisfied in the region 
(i/ — l)/(i/ + 1) < V2/V1 < 1. To choose the appropriate branch we find the 
adiabatic speed of sound in the region behind the shock, (0^)2 ( it is obvious 
that in the region 1 of the cold gas the speed of sound is 0). 

From (jllj) follows 

(c.)^ = {v-m- (|^).] = {v- 1)(1 - n2V2) (15) 

where we use (fTUIl. We express this speed of sound in terms of V2. From ^ 
and 

V 1 

w = -P2 + TF 

z/ — 1 V2 



which follows from (jlip ] we find 



V2 = W2vi = -^j^P2Vi + V2. 

Solving this equation with respect to V2 we obtain 



Therefore 

... V2 2z/ P2V2 



Inserting this in (fT^ and using (fT^ we get 



= (17) 

Now it becomes clear that for reahstic values of > 1 

— < (z/-l) 

(that is (cs)2 > 0) chooses the physical branch of the adiabatic (IHj), that is 
only the left part of the Fig.l has a physical meaning. 

We also demonstrate that the value of parameter v cannot exceed 2, since 
otherwise the flow becomes non-causal. Actually, from the boundaries 

v+l- V^-^ ' 

(where the left inequality has been found above) of the physical region and 
relation (fTTjl follows that the causality condition {cs)\ < 1 imposes the fol- 
lowing restrictions on u 

{u - if <'^^-^{2u -I), that is l<u<2 (18) 

2.3 RATIO OF VELOCITIES BEHIND AND IN FRONT 
OF THE SHOCK 

Here we derive equation for x = which is derived from the condition 

of continuity of the particle flux density 

U2_V2 

Using ()16|) and the adiabatic (fT^ we obtain for its physical branch the fol- 
lowing expression: 

Uo (z/ - 1) - ^ 

_i . ^_L^ (19) 
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Figure 2: Ratio x = P2/ Pi as a function of for 5 values of the four- velocity Ui 
from .1 to 3 



We also provide the expression for x = P2/P1 found in [H],[II]: 







(20) 



Its imphcit solution is: 



u = l + Ulx 



x-1 



(21) 



ul{x-i) + ^i + ul{i-x^) 



The respective graph of x = /(z/, Ul) is shown in Fig.2. 

As is seen from the figure, for 5 values of z/ < 2 ( which ensure the causality 

of the flow) parameter x is always less than u — 1. 

For two asymptotic cases of 

i) a non-relativistic flow {Ui ^ Pf << 1) and 

ii) the ultra-relativistic flow ( with Ui >> 1) 
expression (j2ip is signiflcantly simplifled. 

i)Ui ^ Pi << l,e ^ n, V2 ^ V2. This allows us to rewrite (fT^ as follows 



U2 



x[i-^{i-x')] + o{Pt) 



z/ - 1 - X 



(22) 
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By representing 

X = Xo + XiPl 

we obtain from (j22p the final expression for x 

Using (j2ni) and the fact that V2/V1 V2/V2 = x iii the shock adiabatic (fT^ 
we obtain: 

2 2 
P2V1 = — —v^ 
y + 1 

that is the result obtained in newtonian gasdynamics,(e.g., PP). for a strong 
shock wave. This means that the lower part of the shock adiabatic ( close 
to V2/V1 = (z/ — l)/(z/ + 1), which seemingly contradicts the continuity of 
the particle density flux, actually represents the non-relativistic limit of the 
strong shock wave in Bethe- Johnson gas. 

u)Ui » 1). Solving Eq.^we obtain 



which means that the velocity behind the shock 

(32^iy-l 

On the other hand, 

Using in we find with the same accuracy 



| = X-=(.-l)(l-j^X^) + 0(l/C/f) (26) 



Upon substitution of p6|l in shock adiabatic ()14|1 we arrive at the following 
expression: 

712 ^ 2^-/32/(z^-l) 

and 

^ = ^ = (2 - v)Ul > 1 (27) 
10 



It follows then that in this case 



62 > ei, 

which in particular means that (j27p is true even the particle number density 
is not conserved. The latter is easily verified by direct substitution of equa- 
tions of state (|TT|) and (fT^ in the laws of energy and momentum conservation. 

It was suggested in [T2] that for a compression ratio n2/ni = 5 — 20 the 
baryon matter can be transformed into the quark gas. By using these values 
of n2/ni in the shock adiabatic (jl4p we can find the range of the respective 
transition pressures p2- By substituting 

V2 U 111 111 

T7 = H , 

Vi U — 1 712 ^2 

which follows from the definition of V2 and the equation of state (fTT|) . we 
rewrite (fT^ function of (^2/7^2, ly) : 

u^z^ + z^x[l - x{u - 1)] + 2(z/ - l)x^{[l - x{u - l)f - 1 - 

+x^{y-lf\{v-l)x-u-l\ = Q (28) 

where z = P2^i and x = n2/ni. If we fix the value of ei = rii as 133MeV/ fm? 
( based on the assumption that the ratio of the bag constant B entering the 
bag model equation of state and ei is 0.45 jIU] and B ^ 60Mev/ fm^, then 
numerical solutions of (j28|) in terms of p2 in lO^^dn/cm'^ for the following 
values of the parameters 

X = 5, 8, 10; u = 1.7, 1.8, 1.9 

are given in Table 1 



As is seen from the table, the transition pressure for n2/ni = 8 and 10 are 
of the same order of magnitude as the ones found in Ref. [THj . 

Now we consider the shock wave in an observer's frame of reference such that 
the shock velocity is —D, the velocities in front and the back of the shock 
are f3i and (32 respectively. In the shock frame of references these velocities 
P[ 2 become 
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Table I. Transition pressure p2 for various values x = n2/ni and u 



X =n2/ni 


5 


8 


10 


V = 1.7 


3.51 


42.37 


37.88 


u = 1.8 


5.68 


27.31 


49.30 


u = 1.9 


10.67 


34.76 


42.66 



From the conservation of energy and momentum across the shock follows 
that 

m = (30) 

62 - ei 

Using equations of state for Bethe- Johnson gas pijl we obtain that for a 
strong shock {Ui >> 1): 



= {u-l){l- ^y^) + 0{Ur') ^y-l (31) 

Using in we get a quadratic equation with respect to D. The 
physical solution to this equation in terms of the velocities /3i,/32 in the 
observers frame of reference is 

^ _ v/(2 - vfiP, + ^2? + 4[(z/ - 1) - [1 - P2Pi{v - 1)] 



2[l-A/32(z/-l)] 

(2-z/)(/32 + /3i) 



(32) 



If z/ — > 2 then according to (jHH) D ^ 1. The respective graph oiD = f{Pi,P2) 
for u = 1.3 is shown in Fig. 3. In what follows we need to have the ratio 
W2{P2) /W2{0) where W2{P2) is the value of 1^2 behind the shock (in the 
observer frame) and W^2(0) is the respective value for = 0. Since in a 
strong shock e » ra,this implies 

V 

W2 ^ -p2'T-2 

U — 1 

Therefore we get from (jTTj) 

W2{(32) _ [l-/5^/(z/-l)]^,=o 



W2{0) [l-(3'2/{u-l)] 



(33) 
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Figure 3: Shock velocity D in the observer frame as a function of the ratio of 
velocities f3i and (32 in front and behind the shock respectively 



Since for a very strong shock wave /?i — 1 the velocity I3[ also tends to 1 
according to Therefore from (p?T|l follows that — z/ — 1. As a result, 
fl33|) becomes the indeterminacy of the type 0/0. Using V Hopitale rule we 
find (cf.USI) 



where 



9(^2 /df3i 



(34) 



Performing elementary calculations we find 

MP2), 1-/32 



W2iO)''"'^ I + P2 



(35) 



3 ONE-DIMENSIONAL SHOCK WAVE PROP- 
AGATION 

We consider the following problem jS]: 2 shock waves propagate in opposite 
direction in a "gas cylinder" as shown in Fig. 4. 
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Figure 4: Geometry of 1-dimensional shock waves propagation in the observer 
frame of reference 



The distance travelled by the left-bound shock is denoted by d and the re- 
spective distance travelled by the right-bound shock is denoted by /. A gas 
in front of both shocks is assumed to be "cold", that is its pressure is p = 0. 
A compressed gas is taken to be Bethe- Johnson gas, obeying the equation of 
state (jllll . Initially we choose the frame of reference where the gas between 
2 shocks is at rest. 

When the shock wave propagating to the left with a velocity D reaches the 
left end of the cylinder, there begins an outflow of matter. As a result, a 
rarefaction wave has been formed which starts to move to the right with the 
speed of sound. At this time the shock wave moving to the right with a 
velocity D has not yet reached the right end. Therefore the two possibilities 
can materialize: 

a) the rarefaction wave can catch up and overtake the right-moving shock 
wave before it would reach the right end the right end of the cylinder, 

or 

b) the rarefaction wave will lag behind the right-moving shock wave. 

If we use the shock frame of reference ( of the right-bound shock), then the 
gas velocity behind it is D, and in front of it is given by (j29|) . We consider a 
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strong shock, that is /?i ^ 1 which imphes that the respective velocity in the 
shock frame of reference is also f3[ ^ 1. Therefore combining the relation 
which is obtained by dividing energy and momentum fluxes across the shock 
( see [13) 

^ P2 + D^e2 

and the equation of state (fTT|) we obtain the quadratic equation for D whose 
solution ( excluding the trivial root D = 1) is 

D={u-l) /^/"^ ^ (36) 

Upon substitution of expression ()27|) for 

^^/^^^^^-i ^-/I/"(^'-l) >>^ 

( valid for a very strong shock with /3i — >■ 1 and j32 u — l) in (|35|) we obtain 
that 

D 

This means that we can find the minimum length Imin of the "tunnel" such 
that the rarefaction wave overtakes the shock wave jT^I 

where the speed of sound in gas behind the shock {Ps)2 is given by the 
following expression [0] 

{Ps)2 = 0^ - 1 7 7T— T- (38) 

(z/ - 1) + iyp2 n2 



Taking into account that P2/n2 ^ 1, we see that {l3s)2 — ^ ^/u — 1. Therefore 
from equation ^7\i we get 



Imin _ Vl' — I + I 

~J~ ~ 1 - y^T^ 



(39) 



For the values of u = 1.3 — 1.9 the respective Imin/d varies from 3.42 to 
37.97. In general Imin/di is a monotonically increasing function which tends 
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to infinity when z/ — > 2. As was shown in Refs.jSl and ^2] the entropy change 
in the case of / < Imin is 

As ~ + 1 
a 

The problem becomes more involved if l/d > Imin, because the rarefaction 
wave can overtake the shock wave but cannot pass through it. Therefore the 
former is reflected from the latter thus creating the flow region bounded on 
one side by the shock wave moving to the right and on the other side by the 
reflected wave moving to the left. The resulting flow can be described with 
the help of the equation for a general one- dimensional motion of a relativistic 
gas [T] 

where 

y = Ln{W/Wo), 

Wq is the value of W behind the shock when the velocity (32 = 0, 

rj = tanh-^p (41) 

and the generalized potential x is connected to the spatial x and temporal t 
coordinates as follows 

_ dx dx 
X = e ^{—sinhrj — —coshrf) 
oy or] 

— 

t = e~^(—coshr] — —sinhr]) (42) 
oy or] 

By using (|HHj) in pUj) we obtain 



(.-l)(l-f_!)f| + (2-.+ ^e-)§^-^ (43) 

Wq dy^ Wq dy dr]^ 



(.-l)^ + (2-.)-i--i = (44) 



For a strong shock wave Pi ^ — ^, Wq S> l,y >> 1, and equation 
becomes 

d'^X _^ (2 - u)— - — 

The respective boundary conditions are based on the fact that in the strong 
shock wave ( in the shock frame of reference) the velocity of gas behind the 
shock is = z/ — 1 (cf.pB]). Therefore from equation (j29|) follows that the 
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shock velocity D in terms of ( in the shock frame) and (32 = —tanht] (in 
the observer frame) is 

_ dx u — 1 + tanhr] 

~ dt ~ 1 + {u - l)tanhr] ^ ' 

On the other hand, from and 1)411) follows that at the shock 

y = 27] (46) 



Because equation ()27|) is true for very strong shock waves, even if the particle 
number is not conserved, relation ()m|) remains valid in this case. If we use as 
model an ultra-relativistic gas then equation (j2Z|) is exact, although it refers 
to 62 instead of In addition, the variable y becomes 

T 

y = Ln{ — ) 

where T is the absolute temperature of the gas. Equation is replaced in 
this case by the following 

u — 1 , . 

y = 2^^^ (47) 

where u is the arbitrary adiabatic index. 

If in our problem (considering Bethe- Johnson gas) we assume that the par- 
ticle numbers are not conserved then for a very strong compressive shock ( 
W ^ p/n » 1) thermodynamic calculations show that now 

T 

y ~ Ln{—) 

and 

- - i^y^^'^'^ (48) 

where a is the entropy per unit of volume and we use p We can combine 
both expressions and for y into one generalized relation 

y = 5r] (49) 
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where S can be either greater or less than 1. 

We substitute (jl^ in use (jl^ and equation (jH)) and obtain ( after 

rather lengthy calculations) that at the shock ( where y = 5rj) 

("^ + ^)| + (^ + <l(l-|w..., = (50) 

The other boundary condition (at the Riemann wave) is borrowed from 
Ref.0 

By introducing new variables 



5' x/F^T) 
we rewrite equation (jUJand the boundary conditions (j^Ujl . fj^ 



(5 + 1 ^ + ;/- 1 g b'^ -nu^ \ d \ d 1 g ^ 

(53) 

Applying Laplace transform to variable h in equations ()51 |) -()53 |) we obtain 



( J ;a ? ~ IJtt + 2(1 J )r —\— + —==rx = (54) 

and 

,,6+1 6 + 1^ — 1, 6'^ ~ nu + 1 d Id 1 

- 1 - 1 d{u - 1) dq dq ^v-\ 

where r is the transformation parameter. 

We represent a solution to (jSD), as a{r')e'^^^~^^ . This leads to 2 algebraic 
equations: 

(1 - 52)P + A;[2(5 - l)r - + = (56) 
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^ -f)il + k + f)a+ ^ ^ A = (57) 

5+1 ^/v-l 6 + 1 

where ^ 5 \ 

k = f = , 6 = , X 



y 7^7^' 

In turn, the variable A is calculated at the point where Riemann's wave 
catches up with the shock wave |3] 

dx dx 
^ ^ ("^)l<?=o,fc=o = -(^)|g=o,/i=o 



We will be interested in finding the entropy change in this fiow. It is given 
by the following expression jlj 

= Eo / au2dt (58) 
Jti 

Here Sq is the "tunnel's" cross-sectional area, a is the entropy density behind 
the shock, U2 = (^2/ a/I — ti is the moment of time when Riemann's wave 
catches up with the shock wave and ^2 is the moment of time when the shock 
wave reaches the right end of the "tunnel" (Fig. 5). Since n we consider 
a strong shock P2 ^ ^ 1, that is U2 ^ — l)/-\/2 — i'- Because the 
time interval dt is measured in the shock frame reference the respective time 
interval in the observer frame ( where initially both waves, the strong shock 
and Riemann's wave move with a velocity D) is dt/ \/l — D"^. Taking these 
facts into account and using and PH|) we obtain from (|HH| 

where Uc is the value of y at the moment when the shock wave reaches the 
right end of the " tunnel" . 

By introducing variables q and h into (j59p and effecting the integration we 
find after rather involved calculations 

AS = ai([l - /i'e"2 ± ase^'^ / ' \ ~ Z±{h'^ - z)dz] (60) 

Jo '^c~ ^ 



19 



where 



^ ±SoaoA {A + + yA^] 

2 A + ^ 1 + 



(A - V^2^1)(i/ - 2^/^7^ 
as = 



2 - z/ 

57(1 - for 5 < 1 
57(5^ - 1) for 5 > 1, 

h 2-1/ 



j/A + 2V(^2^1)(^^-i: 



2 — 1/ 

/i(/i'^ — z) Bessel function of the imaginary argument 
Ja{h'c — z) Bessel function of the first kind, 

and everywhere the upper and lower signs correspond to A > 1 and A < 1 
respectively. 

We also have to find the value of he entering integral in (jUUI) . To this end we 
use the relation from Ref. [H] 

Xc-l = to-tc (61) 

where to = A/a/z/ — 1 is the moment of time when Riemann's wave overtake 
the shock wave ( in the observer's frame), L is the coordinate corresponding 
to this location of the wave, Xc and tc are the coordinate and moment of 
time describing the state of the shock weave at the right end of the "tunnel" . 
Inserting x and t from (jl^ in (jUT|) and performing rather lengthy calculations, 
we obtain 

+ \){y/17^ - u + 1) _ 



^^'^ ^^'' {A + VW^exp[-2uAh'/{2 - u)] + A - V^^l} - 



2A 



/ Z^{K-z)dz (62) 
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where 

z/A-2v/(A2^1)(z/-l) 
2 — u 

In what follows we neglect the integral in (jU^ which would introduce an 
error not exceeding 2%. 

From the kinematic considerations follow that 

A ,, I - d 
L + , = {l + d) + 

and 

Therefore relation ()62|) (with the integral dropped) would yield the expression 
for exp{—2Ah'/i') in terms of u and A 

±iA + v/I^)2exp(-^^) = 
[uA{l - 7^7^)3 + ^ l)(z/ - 1) - A(z/ - V^7^)]{63) 



^fV^[^A - ^ 1) ' ci 



i) If (5 = 2(z/ — ( as in the case of an ultra-relativistic gas with an 
arbitrary constant adiabatic index v) implying A = — 1/ (2 — v) then 
we choose in (p]j) the lower sign. ^ As a result, we obtain from ( IUU|) with 
the help of dHl 

^ 3z/-2 + (z/ + 2)v/^7^ 2-z/ , 2z/ / _ ^..(3.-2)72. _ 

SoAao 2(3i/-2) + 27^7^^ ^2 - z/ ci ^ 

(64) 

where 

l_ ^ 1 + y/^T^ 

d - \ - 

The respective graph of < AS >= AS*/ (SqActo) is shown in Fig.5. Expression 

^We could have chosen the upper sign. However, the requirement A > \ would have 
left us with a very restricted range of parameter v.l <v< 4(1 — ■\/l/2. 
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Figure 5: Variation of the dimensionless entropy change as a function of the 
adiabatic index v and Ijd 



given by yields as a particular case the expression found in |^ for z/ = 
4/3 It is seen that the entropy change in region 2 increases and reaches its 
maximum for the values of v in the range from 1.3 to 1.4. A further increase 
of V results in a sharp decline of the entropy change, approaching the purely 
adiabatic regime for the threshold value oiv = 2. 

a) If, we consider 5 = 2 ( as in 0^1) we obtain in the same fashion as in z: 
^ 3v/z/- l + z/+l _ 2u l/d _ ^.p(^_i)/2^. 

(65) 

where Im/d is given by (p^ . From follows that in this case the entropy 
change can be positive only if u > 2. However this also means that the flow 
becomes non-causal , since ( as we showed above) in this case jSs > 1. There- 
fore if 5 = 2 the one-dimensional flow of Bethe- Johnson gas with u < 2 is 
characterized by a zero change of entropy when Riemann's wave catches up 
with the shock wave (and stays zero from that moment on), that is the flow 
becomes adiabatic. This contrasts with the previous case i) which is a gen- 
eralization of the analogous flow for an ultra-relativistic with an arbitrary 
adiabatic index ( not necessarily i/ = 4/3 as in ^). In the latter the entropy 
change is linear in l/d before Riemann wave catches up with the shock wave, 
and obeys the power law for l/d (given by eq l65|) after Riemann wave over- 
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takes the shock wave. 



The flow described here can be associated with a nucleon-nucleon colhsion 
(cf.y^) which is modelled as a collision of the nucleon with a cylinder cut 
out of the nucleus and whose cross- sectional area is equal to the one of the 
nucleus. The nucleus linear size varies from the diameter of the nucleus to 
the "diameter" of the nucleon, d. It is clear that such an approximation 
can be used only to obtain certain estimates about parameters of the real 
process, for example to evaluate an entropy change. The ratio l/dis related 
to the atomic number A: l/d = A^/^. If we take max A = 257 this would 
give us the maximum value of z/ = 1.53 ( according to eg ]39|) which is close 
to the region where entropy change reaches its maximum. 

4 CONCLUSION 

We found the shock adiabatic for BJ gas and applied the found detailed 
relations for a study of a one- dimensional flow of such a gas in a context of 
a nucleon-nucleon collision. If baryon matter can be transformed into quark 
matter at a certain range of the compression ratio (5 — 10), then the proposed 
model allows one to identify the baryon matter as Bethe- Johnson gas with 
the adiabatic index u ^ 1.5. 

Such an approach generalizes an analogous approach based upon an ultra- 
relativistic gas. The latter represents a particular case of BJ gas with adia- 
batic index u = 4/3. The entropy change behind the shock strongly depends 
on the choice of the adiabatic index in BJ, by reaching its maximum in the 
range of values of u from 1.3 to 1.4. 
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